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Abstract: We discuss in detail the possibility of observing pairs of simultaneous 
parallel air showers produced by the fragments of cosmic ray nuclei which disinte- 
grated in collisions with solar photons. We consider scenarios with different cosmic 
ray compositions, exploring the predicted rates for existing and planned detectors 
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j§ \ and looking for methods to extract information on the initial composition from the 

characteristics of the signal. In particular, we find that fluorescence detectors, such 
as HiRes or the Telescope Array, due to their low threshold (~ 10 17 eV) and large 
area (~ 10 4 km 2 ) may observe several events per year if cosmic rays at those energies 
are indeed heavy nuclei. The possibility of exploiting the angular orientation of the 
plane containing the two showers to further constrain the cosmic ray composition is 
also discussed. 
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1. Introduction 

A long time ago, Zatsepin and Gerasimova [p], suggested that the photo-disin- 
tegration of cosmic ray (CR) nuclei with solar photons could allow to study the 
CR composition at very high energies, E > 0.1 EeV (1 EeV= 10 18 eV). Indeed, if 
a nucleon is stripped from a heavy nucleus on a photo-disintegration process, the 
comparison of the energies of the surviving nucleus with that of the emitted nucleon 
would directly provide the mass of the primary nucleus. 

In order to be able to observe this kind of events, one needs that the probability 
for the disintegration process to take place as the nucleus traverses the solar system 
be non-negligible. It is also necessary that the separation among the fragments be 
sufficiently large so as to allow their individual identification but however not larger 
than the spatial coverage of the CR detectors. 

To maximise the photo-disintegration probability, previous works have focused 
on CR energies for which the interaction with solar photons takes place at the giant 
resonance for photo-disintegration, i.e. such that the typical photon energy in the 
CR rest frame is 10-30 MeV (e.g. E ~ 1 EeV for iron nuclei). 

Regarding the separation among the fragments, Gerasimova and Zatsepin (GZ) 
originally estimated that it was of the order of 1 km 0], based on the splitting 
resulting from the transverse momenta acquired by the fragments in the disintegra- 
tion process. However, as Zatsepin later realised, the dominant contribution to the 
splitting actually results from the deflection of the fragments in the solar system 
magnetic field (which was initially largely underestimated). This deflection depends 
on the charge to mass ratio of the fragments, and is hence different for nuclei, protons 
or neutrons. 



In a recent reanalysis of this mechanism, Medina- Tanco and Watson J| have 
performed detailed orbit integrations for incident Fe nuclei, using a realistic model 
for the solar system magnetic field. With these improved computations, they showed 
that the resulting separations are typically much larger than 10 km, and hence exceed 
the typical sizes of existing detectors, what strongly limits their observability. Going 
to higher energies in order to reduce the average shower separation down to an 
observable range has however the problem of a much reduced rate, since not only the 
cross section decreases outside the giant resonance, but more importantly the CR 
flux falls abruptly. Hence, the conclusion of that work regarding the observability of 
the effect was not encouraging. 

In this work we want to extend previous analyses in several ways. Since the 
ultra-high energy CR are most probably not dominantly Fe nuclei, we study the GZ 
effect for different initial compositions. In order to cover all the range of possible 
CR nuclei masses we take as prototypical examples of heavy, intermediate and light 
compositions the nuclei 56 Fe, 16 and 4 He. We find that although the cross sections 
are reduced going to lighter nuclei, for a given initial energy the separation among 
the fragments are smaller (approximately proportional to the mass number A), and 
hence this helps to make this kind of events observable. Also, since the energy of the 
emitted nucleon, which is smaller than the parent cosmic ray energy by a factor A, 
becomes higher the lighter is the initial mass composition, this allows for an easier 
detection of the nucleonic fragments. 1 We have also extended the cross sections 
beyond the pion-production threshold in this study. 

The determination of the parent CR mass through the comparison of the ener- 
gies of the two showers may be problematic in some cases. This occurs for instance 
if the determination of the nucleonic fragment energy is poor, as would happen if 
the energy is small so that only one detector in a ground array is hit by the shower 
(e.g. for a fragment energy ~ 10 17 eV for Auger), or if more than one nucleon is 
emitted in the photo-disintegration. One may then attempt to use the size of the 
separation among the fragments to infer the initial composition. However, for a 
given initial energy the separation among the fragments also depends on the dis- 
tance from the Earth at which the photo-disintegration took place. Hence, nuclei 
with different masses but interacting at different distances may give rise to equally 
separated showers. We have found that this degeneracy can be lifted in many cases 
by looking at the inclination of the plane containing the two fragments, since this 
angle depends on the distance from the Earth at which the photo-disintegration 
took place. Looking then at the whole available information (total energy, separa- 
tion and inclination of the fragments) it will be possible to get further insights into 
the initial CR composition. 

1 For a Fe nucleus with E ~ 7x 10 17 eV, as considered in jjj, the emitted nucleon has E ~ 10 16 eV, 
which is below the threshold of the large extended air shower arrays. 



2. The photo-disintegration process 



As a CR traverses the solar system, it encounters a flux of energetic (e 7 ~ 1 eV) 
solar radiation. The lifetime for photo-disintegration off this radiation, computed in 
the rest frame of the CR nucleus (primed quantities), is just 

I roc Hn' 

- = c/ o de<^( £ <), (2.1) 

with er(e') the photo-disintegration cross section. Transforming back to the lab 
(Earth) frame, using that e 1 = £7(1 + j3cosa), where (5 — v/c ~ 1 in terms of 
the CR velocity v, 7 is the usual relativistic factor and a is the angle between the 
directions of propagation of the CR and of the photon, one finds for the CR mean 
free path 

W) = r de te i0a{e ' )[i+(3c ° sai (2 ' 2) 

Here £ is a coordinate measuring the distance from the Earth along the arrival 
direction of the CR, so that if £ is the unit vector in that direction, the CR velocity 
is v ~ — c£ and cos a = £ • f, with r the spherical radial coordinate centered in 
the sun. 

The probability that an incoming CR suffers fragmentation along its path to- 
wards the Earth is then 
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This probability turns out to be small (<C 1CT 3 ), since A(£) is typically much larger 
than the characteristic solar system dimensions. For practical purposes, we will 
integrate up to £ maa; = 5 AU, since the contribution from larger distances turns out 
to be negligible due to the decreasing photon flux. 

The solar photon flux can be obtained as a blackbody spectrum with the tem- 
perature of the solar surface, T s = 5770°K (/c_bT s ~ 0.5 eV), normalised so as to 
reproduce the solar luminosity, L & = 4nr 2 c j de edn/de. One then finds 

^ = 7.2X10' — (l^V [eV cm]-' . (2.4) 
de exp(e/A; B T s ) - 1 V r J [ J V ; 

Regarding the photo-disintegration cross section a, one can distinguish different 
regimes according to the photon energy in the nucleus rest frame. There is first the 
domain of the giant resonance, from the disintegration threshold (e' ~ 2 MeV) up 
to ~ 30 MeV, in which a collective nuclear mode is excited with the subsequent 
emission of one (or possibly two) nucleons. Beyond the giant resonance and up to 
the pion production threshold (ef ~ 150 MeV) the excited nucleus decays dominantly 
by two nucleon (quasi-deuteron effect) and multinucleon emission. Detailed fits to 



the cross sections and branching fractions in these first two regimes were performed 
by Puget, Stecker and Bredekamp j| in their study of the CR photo-disintegrations. 
For our purposes, it will be enough to adopt the simpler expressions for the total 
cross section given in [|], which are 

<K0 = *g*(0 = 1A5A TT2 ir77™ mb fOT e ' < 30 MeV > (2-5) 

[e l - eg) 2 + [el y 

where T = 8 MeV and e = 42.65A" - 21 MeV for A > 4 and e = 0.925A 2 - 433 MeV 
for A < 4. For higher energies we take 

o-(e') = max{a GR (e'), A/8 mb} for 30 MeV < e < 150 MeV . (2.6) 

We do not cutoff ctcr above 30 MeV since for light nuclei the peak of the giant 
resonance is already close to 30 MeV and the resonance is wide. 

When considering light nuclei and/or very energetic CR (E > 1 EeV), it is 
necessary to consider in some detail the cross section above the pion production 
threshold. From this threshold up to e' ~ 600 MeV, the disintegration is dominated 
by the delta production ( / ~fN —>■ A —>■ Ntt, where N is a nucleon). 2 In this regime, the 
cross section per nucleon is found to be almost independent of the target nucleus || 0]. 
Due to nuclear Fermi motion and interaction effects, the peak in this cross section 
is reduced and the resonance becomes somewhat wider than in the free proton case. 
(Higher resonances beyond A(1232) should also contribute for e' ~ 1 GeV.) 

We have parametrised the total photo-absorption cross section in the region 
above the pion production threshold as 



a(e') = A 



1 Q ~ /I 
- + b e exp 



mb, for e' > 150 MeV , (2.7) 



where e = (e' — 150 MeV)/ei. A good fit to the experimental data |J is obtained 
with ei = 180 MeV, S = 0.3 and v = 1.8, so that the maximum in the cross section 
(o 'max I 'A = (l/8 + S) mb~ 0.4 mb) occurs for an energy e' = ei + 150 MeV~ 330 MeV. 
The value of v gives the right size of the resonance width. 

Near the resonance peak, the one nucleon emission by direct knock-out is rele- 
vant, but the nucleon multiplicity increases at higher energies. The pions emitted 
will decay, for the typical 7 factors considered (7 ~ 10 7 -10 9 ), in a distance much 
shorter than 1 AU, and only their decay products can reach the atmosphere. How- 
ever, these decay products will have an energy much below the initial CR energy and 
hence will pass largely unnoticed. Furthermore, the muons from -n^ decays will be 
significantly deflected by the solar system magnetic fields and hence arrive far from 
the nucleonic fragments. 

2 We note that using Rudstam type fits to the proton-nucleus cross section [|[ [)) to estimate the 
photo- nuclear rates, or also photo-disintegration data from bremsstrahlung photons |To| |, one misses 
the strength of the A resonance peak. 



3. The deflection of the fragments 



After the parent CR photo-disintegrates, the fragments travel towards the Earth 
and are deflected by the solar system magnetic field with respect to their initial 
trajectories. The resulting deflection can be computed by integrating the transverse 
displacement due to the Lorentz force from the fragmentation point, at a distance £, 
up to the Earth. 3 In this way one finds 

d '«> = ^M>' B «"> x «- < 3d > 

where e and m are the proton charge and mass. The deflection df of each fragment is 
then proportional to its charge to mass ratio, Zf/Af. For the separation among the 
fragments what matters is the difference \Z\j A\ — Z 2 /A 2 \. If the emitted nucleon is 
a neutron, this is just Zj (A — 1), with Z and A associated to the charge and mass of 
the parent nucleus. If a proton is emitted, the difference is given by (A — Z) / (A — 1) . 
Thus, the separation among the two fragments is approximately given by eq. ( |3.1| ) 
but replacing Zf/Af by 1/2 for Fe and O, and by 2/3 for He. In the case of multiple 
nucleon emission, all protons will be deflected by the same amount by the magnetic 
field (and all neutrons will be undeflected), so that they will produce essentially 
one shower (the typical separation of the strongly overlapped nucleonic showers, 
see footnote |3], is smaller than the separation of the detectors in the largest arrays, 
which is e.g. 1.5 km for Auger). The interesting possibility also appears of having 
three simultaneous showers caused by the main fragment, the emitted proton(s) and 
neutron(s), all contained in the same plane and equally separated. 

As in ||, we consider for the solar system magnetic field the realistic model of 
Akasofu, Gray and Lee [IT], which describes the field within the heliosphere, taken 
as a sphere of radius r 2 = 20 AU. It consists of four components (written here for an 
odd solar cycle, i.e. with the north polar region having the S magnetic pole): 

a) The solar dipole contribution can be written, in cylindrical coordinates around 
the sun, as 

BfP = -d(p 2 - 2z 2 ){z 2 + p 2 )- 5 / 2 , (3.2) 

< P = -^P^ 2 + pT 5/2 , (3.3) 
Bf p = , (3.4) 

where B s r\/2 is the magnetic dipole moment of the sun. Taking ri to be the 
solar radius, r\ = R & , one has B s = 2 Gauss. 



3 In the following we neglect the deflection due to the transverse momenta acquired by 
the nucleonic fragments in the photo-disintegration, since it leads to separations < 0.2 km 
y/A/(E/ EcV)(£/l AU). The corresponding separation for an emitted alpha particle would be 
four times smaller. 



b) The dynamo component arises from the rotation of the sun in the dipole field. 
It is generated by a sheet current distribution that flows outwards along the 
polar axis, then flows along the heliosphere towards the ecliptic plane and then 
radially inwards back to the sun. The resulting field is in the <fi direction, 



Bf n = sign(z)B 4 
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where B^ = 3.5 x 10 5 Gauss and po = 1 AU. 



c) The ring current component arises from a sheet equatorial current extending 
up to 7*2, with current density oc p -2 . The resulting magnetic field is given by 
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Jo,i are Bessel functions and B po = —3.5 x 10 5 Gauss. 

For the region of the solar system interesting for the present work (?*i r < 
5 AU) the integrals in eqs. ( |3.6| ) and ( |3.7| ) can be approximated by taking the 
limit t*i — ► and r 2 — > oo in G(k). This gives 
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(3.9) 
(3.10) 



d) The sunspot component is approximated by an ensemble of 180 dipoles of 
radius 0.1 R Q located in the equatorial plane, just below the solar surface at 
a radius of 0.8 R Q (their expressions are similar to eq. (|3.3j ) with B d = 1000 
Gauss). They provide a strong magnetic field near the sun surface and allow 
the ring current field lines to connect to the solar surface. 

All the components change sign with the eleven years solar cycle. We plot in 
fig. [I] the three components of the magnetic field B^, B p and B z as a function of 
p for three different values of z. At large distances the dominant contribution is 
given by B^ and B p (due essentially to the dynamo and ring current components 
respectively). 4 

4 The B z component is at variance with the one plotted in ref. and is proportional to r~ 3 for 
large radius. However, being this one the smallest componet, it has little impact on the results. 



6 



o 
o 

o 

CO 



10000 
1000 
100 
10 

1 

10000 
1000 
100 
10 

1 

10000 
1000 
100 
10 

1 




z = 0.1 AU 




----_+ 



z=1 AU 



+ 



r[AU] 

Figure 1: The three components of the solar system magnetic field, Bs (solid line), B p 
(dotted line) and B z (dashed line), as a function of the radial cylindrical coordinate p for 
values of the z coordinate above the ecliptic z = 0.01 AU, 0.1 AU and 1 AU. We have 
plotted the modulus of each component and the sign is indicated in each curve. 



We notice that the field model adopted gives a consistent description of Parker's 
spiral magnetic field configuration, and allows us to discuss the general features of 
the GZ events. Some departures from this simple model occur for instance due to the 
variation of the overall field normalization during the solar cycle or from a possible 
tilt of the magnetic equator with respect to the ecliptic, and these can be taken into 
account in the event of an actual detection. Sporadic magnetic perturbations such 
as shock waves produced by solar flares could also add some noise. 

We describe the CR position by the distance £ to the Earth and the latitude and 
longitude angular coordinates (b,£), defined analogously to the Galactic ones, with 
b = in the ecliptic plane and the sun at the origin of £, with increasing £ values 
to the left (with the N up). The separation among the fragment positions can be 
deco mposed along the b and £ directions, d = dgi + dbb. We will use the modulus 
d=Jdj + d\ and the angle with respect to the £ direction, atan(d 6 /(^), to describe it. 



We can use the symmetry properties of the magnetic field about the ecliptic 
plane to infer the relationship among the separations between the fragments arriving 
from directions with the same longitude £, but opposite latitude b, fragmenting at 
the same distance from the Earth (so that the sign of £ 2 changes, but the component 
of £ in the ecliptic plane is the same). Under the reflection z — > —z, the radial 
and azimuthal components of B change sign, while B z is invariant. Thus, the z 
component of the product B x £ will change sign, while the orthogonal component 
remains unchanged. This means that db(—b,£) = —db(b,£) and di(—b,£) = de(b,£). 
Or, in other words, rays arriving from the same £ but opposite b and produced at the 
same distance from the Earth have separations between fragments equal in modulus 
but with opposite angles with respect to £. Hence, when integrating the fragment 
orbits, only half of the celestial sphere needs to be evaluated. 5 We also note that in 
an even solar cycle all the field components are reversed, but this leaves unchanged 
both the separations and the angles. 

4. Results 

Typical separations for different parent nuclei are illustrated in fig. |^ for two dif- 
ferent energies. The mean separation for each arrival direction is plotted (dashed 
lines), where the mean is weighted with the fragmentation probability for the photo- 
disintegration processes to occur at different distances from the Earth, i.e. 



where dr] GZ /d£, ~ 1/A(£). 

For Fe nuclei with E = 1 EeV most of the separations are larger than 100 km 
and are hence out of the reach of existing detectors, as noticed in ref. J|. For 
lower energies this problem worsens, and for instance for E = 0.1 EeV most of 
the separations are larger than 500 km for Fe. The situation improves for lighter 
nuclei, due to the factor 7 in eq. (|3.1|) , which makes the separation approximately 
proportional to A for a fixed energy. For example, for He nuclei with E — 1 EeV most 
of the separations are smaller than 50 km, and hence within the range of existing or 
planned detectors. 

The solid contours in fig. ^ correspond to the fragmentation probabilities rjaz 
(and are labeled by rjcz x 10 6 ). Values of rjcz ~ 10~ 5 -10~ 4 are obtained for Fe nuclei, 
while rjcz ~ 10~ 6 -10~ 5 result for He and O. 

In order to compute the rate of observation of these kind of events in a detector 
covering a surface S, one must only consider showers separated by distances smaller 
than d max ~ y/S. The rates will be a convolution of the CR flux, the fragmentation 

5 This symmetry is not apparent in ref. [3], probably due to insufficient numerical accuracy. 
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Figure 2: Mean fragments separation (dashed curves) and fragmentation probability (solid 
lines) contours for total energy E = 0.1 EeV and E = 1 EeV and for heavy (Fe), interme- 
diate (O) and light (He) nuclei. The separation contours are labelled by their values in km 
units and the probabilities in units of 1CP 6 . 

probability and the fraction of GZ events with d < d max , fd max (E). The areas (and 
thresholds) of the large extended air shower detectors we will consider are 100 km 2 
for AG AS A |2| (E th ~ 1 EeV), 3000 km 2 for Auger @ (E th ~ 1 EeV), fewxlO 4 km 2 
for HiRes |14| or the proposed Japanese Telescope Array |R| fluorescence detectors 
(E t h — 0.1 EeV) and ~ 10 6 km 2 for the proposed OWL |16| type satellite fluorescence 
detectors (E th ^ 10 EeV). 

Regarding the CR flux, above the knee in the spectrum (Ecr > 3 x 10 15 eV), 
the measured flux is 

(4.2) 



$(E C r >E)^A7 (^jp)' [km 2 yr sr]" 1 



Due to the steepness of the spectrum (d<&/dE oc E 3 ), we will just estimate the rates 
of events from CR with energy above a given value E as 



Rate(E CR > E) ~ $(E CR > E) VGZ (E)f dmax (E)Se dc n 



(4.3) 



where fluorescence detectors have typical duty cycles e^c — 0.1 and Q is the solid 
angle. 
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Figure 3: Expected rates of events per year and sr for a detector area of 1000 km 2 and 
for separation between fragments smaller than 1000 km (upper curves), 100 km (middle 
curves) and 10 km (lower curves). Solid lines correspond Fe, dashed lines to O and dot- 
dashed lines to He. Fig. 3a is for a night-time direction (6 = 0°, I = 150°), while fig. 3b 
for a direction close to the sun (6 = 0°, £ = 30°). 



In fig. |3| we plot these rates as a function of the energy of the parent CR, for 
the three nuclei and for maximum distances of 1000, 100 and 10 km. Fig. is 
for a night-time direction (b = 0°, i = 150°) while fig. is for a direction close 
to the sun (b = 0°, £ = 30°). 6 It is clear that, if d max < 100 km, the maximum 
rates per unit coverage are attained in the night side for energies E ~ 10~ 2 -1 EeV 
(depending on the nucleus considered) and exceed one event per year for coverages 

6 Note that in the day-side, \£\ < 90°, for increasing \b\ the average deflections decrease, enhancing 
the rates with respect to the 6 = results of fig. pi For the night-side the opposite will happen. 
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S Q edc 10 3 km 2 sr. Note that the AGASA array has an order of magnitude 
smaller coverage, while the Auger array, with S ~ 3000 km 2 looks in this sense more 
promising. However, to be able to observe these events, low thresholds are required, 
and the emitted nucleon will not fire a tank in Auger if its energy is below 0.1 EeV. 
This means that to see both fragments, the parent CR should have an energy larger 
than 0.1 A EeV, which for incident Fe nuclei is much larger than the typical threshold 
of the detector (Eth — 1 EeV, corresponding to triggering with ~ 5 tanks). Hence, 
the rates for Auger are also quite small (< 10 _1 events/yr). The most promising 
experiments for detecting these kind of events seem to be the large fluorescence 
detectors in construction (HiRes) or planned (Telescope Array), since besides having 
large areas (> 10 4 km 2 ) they will have low thresholds (0.1 EeV). Furthermore, 
it is even possible that if the heavier fragment is above the threshold and triggers 
the detector, the lighter fragment is seen even if it has an energy somewhat below 
0.1 EeV. With these arrays one will then expect few events per year for energies below 
1 EeV if the CR are indeed not mainly protons. Larger detectors, as the proposed 
OWL and AIRWATCH satellites, have much higher thresholds, E th ~ lOEeV, and 
hence will have small rates (< 0.3 events/yr). 

Considering the day side (|£| < 90°), only ground arrays can be employed, and 
looking at fig. |3|& we see that the rates are very small for d max < 100 km. This is 
due to the fact that most fragment separations are larger than the typical size of the 
detectors and become hence unobservable. 

In the case of a positive detection of two simultaneous showers, if a single nucleon 
was emitted, the initial CR composition could be determined as A ~ (E 1 + E 2 )/Ei, 
with E\ the energy of the less energetic shower and E 2 that of the most energetic 
one. This procedure will not lead to a precise determination of A if the nucleon 
shower energy is close to the threshold, and hence poorly determined, or if there 
was multiple nucleon emission, so that the smaller shower is actually produced by 
an unknown number of nucleons. There is however additional information about the 
CR composition, encoded in the separation between the detected fragments, which 
can be used to obtain a good determination of A also in these cases. From fig. it is 
clear that, for a given initial energy and arrival direction, the average separations are 
quite sensitive to the CR composition. However, the distributions of separations are 
quite flat (see fig. ^) due to the wide range of distances from the Earth at which the 
photo-disintegration can occur (with the exception of the directions close to the sun, 
e.g. fig. |]a, for which the most likely location for the photo-disintegration interactions 
are for £ ~ 1 AU). This can result in different nuclei producing similar separations 
if they disintegrate at different distances. 

We have found that this degeneracy can be lifted in many cases by looking at 
the inclination of the showers (i.e. at the angle atan((4/c^)). In fig. |5| we show the 
relations obtained between the separation d and the inclination angle, for b = 10° 
(fig. §a) and for b = 45° (fig. §6) and for different values of £ = 30°, 60°, 90°, 120°, 
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Figure 4: Distribution of the separation between fragments for a total energy of E = 1 EeV 
and for the three nuclei considered, Fe (solid line), O (dashed line) and He (dot-dashed 
line). Fig. 4a is for (6 = 0°, I = 30°), while fig. 46 for (6 = 0°, £ = 120°). The histograms 
give the probability that the separation falls in the range of distances corresponding to 
each bin. 



150°, 180°, 240° and 330°. One moves along the curves as the photo-disintegration 
distance £ is varied. For the southern hemisphere (b < 0°), the inclination angle 
changes sign. Since the separations scale as 7 _1 oc A/E, what is plotted is d x 
(A/56) x (EeV/E). 

Hence, knowing the total energy E, one may infer A by using plots like those in 
fig. This is just done by plotting the d vs. angle curve for the arrival direction 
(b, tj of the observed CR. Confronting then this curve with the measured separation 
and inclination of the showers, and knowing the total energy E, A can be determined. 



12 





Figure 5: Separation vs. inclination angle of the fragments for arrival directions with 
b = 10° (fig. 5a) and b = 45° (fig. 56), and for the £ values indicated in each line. 



One can see that for many directions (e.g. 0° < I < 180° in fig. [|a) the inclination 
of the showers varies considerably as £ is varied, and hence the inclination contains 
the information on the distance to the photo-disintegration point which allows the 
determination of A. In other directions (i.e. those leading to almost vertical curves) 
the resolving power of this method is instead not good. 

In summary, we have explored in detail the possibility of detecting double shower 
events originating from the fragments produced in the photo-disintegration of a CR 
in a collision with a solar photon, for different compositions of the primary CR. 
The best possibilities are for the future large area fluorescence detectors, like HiRes 
and the Telescope Array, for which few events per year are expected if the CR of 
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energies larger than 0.1 EeV are indeed nuclei. These will be rare, peculiar events, 
but it would be interesting to detect them as they can give important information 
about the CR composition at ultra-high energies. The proposed method of using the 
separation and inclination of the showers to complement the analysis with the ratio 
of energies will also help to get a more precise determination of the composition of 
the parent CR. 
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